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We study the stochastic behavior of a set of chaotic vortex loops appeared in imperfect Bose gas. 
Dynamics of Bose-gas is supposed to obey Gross-Pitaevskii equation with additional noise satisfy- 
ing fluctuation-dissipation relation. The corresponding Fokker-Planck equation for probability func- 
tional has solution 'P{{tp{r)}) = 7Vexp(— // {tp{r)} /T), where H {^/)(r)} is the Ginzburg-Landau free 
energy. Considering vortex filaments as topological defects of field tp{r) we derive a Langevin-type 
equation of motion of the line with the correspondingly transformed stirring force. The respective 
Fokker-Planck equation for probability functional 'P({s(^)}) in vortex loop configuration space is 
shown to have a solution in the form of 'P({s(5)}) — J^exp{—H {s} /T), where Af is the normalizing 
factor and H {s} is energy of vortex line configurations. Analyzing this result we discuss possible 
reasons for destruction of the thermodynamic equilibrium and follow the mechanisms of transition 
to the turbulent state 

PACS-number 67.40 

Quantized vortices (QV) play a fundamental role in the properties of superfiuids and other quantum systems. 
Phase transitions, kinetic properties and many other physical phenomena can be effected by QV. For that reason, 
they have been the object of intense study for many years (for review and bibliography see [1]). Studies of the QV 
dynamics have been most successful in simple cases such as vortex array in rotating helium and vortex rings. These 
simple cases, however, are very rare. Due to complex dynamics, initially straight lines or rings can evolve into highly 
entangled structures called vortex tangle (VT), To describe any phenomena involving VT, some information about 
its structure and statistics must be known. One of the most popular ways to overcome that problem is to treat the 
vortex loops as a kind of thermal excitation and to use the thermodynamic methods. The first example of that way 
was the use of the Landau criterium for critical velocity where vortex energy and momentun were used in the relation 
having purely thermodynamic sense. The more recent example is the Kosterlitz-Thouless theory, as well as its 3-D 
variant, currently being developed for review and bibliography see [3]). Besides of vortices in the condensed matter, 
the view of strings as a set of topological defects being in thermodynamic equilibrium, is intensively elaborating now 
in the field theory and in cosmology (see e.g. [4], [5], [6]). 

On the other hand , numerous experiments involving counterflowing Hell, as well as numerical simulation of vortex 
line dynamics [1], [7], demonstrate that the system is not at equilibrium. On the contrary, it possesses features of 
turbulent phenomena. So, for example, the direct numerical simulations ( [8], [9], [10], [11]) demonstrate the presence 
of a cascade transfer of energy both in space of the sizes of vortex loops, and in space of Kelvin waves on every 
individual loop. Furthermore, in the experiments involving vortices in quenched "^He and ^He (see [12], [13]), the 
decay of vortex tangle was observed rather than a thermalization as expected for normal statistical systems. There 
are also serious discrepancies in the statistical description of systems. As follows from theoretical predictions founded 
on the study of dynamics quasi-equilibrium topological defects [5] , decay of the vortex tangle (rate of change of the 
density of vortices, or their total length in unity of volume) occurs by the law 1/t. Though this law agrees with 
the decay of turbulent vortices (see [14]), relaxation of quasi-equilibrium vortices, occurs by some orders faster ( [4], 
[5]). Besides, to describe decay of the turbulent tangle, Vinen has applied arguments used in the theory of classical 
turbulence. Furthermore, distribution of vortex loops with respect of their lengths n{l) = dN(l)/dl in turbulent 
helium obeys to the law n{l) oc (see [15]), while the equilibrium vortices obey to the law n{l) (x (see [16], 

[17]). The very important question is generation of turbulent vortices in the system. In the thermodynamic approach 
vortexes arise due to thermal fluctuations of the order parameter, or remain at prompt quench through the lambda 
line. The mechanisms of appearing of turbulent vortices are not understood. Therefore, the transition phenomena in 
the superfluid turbulence described by the Vinen equation [14], require introduction of phenomenological "generative" 
term whose origin is not clear. 
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Examples listed above show, that chaotic vortex tangle, appearing in quantum fluids can display both the thermo- 
dynamic equilibrium features and the extremely nonequilibrium, turbulent properties. Naturally, there is a question: 
"What state, equilibrium or turbulent, is realized in each concrete experiment, and what sorts of the theory is required 
to use at description of one or another phenomenon?" The answer to this question can give the investigation devoted to 
study of specific mechanisms, responsible for establishment of thermodynamic equilibrium and mechanisms featuring 
destroying the equilibrium and transition to the turbulent state. In this work we mainly concentrate on the first part 
of investigation formulated above and follow how the vortices appeared in a weakly non-ideal Bose gas are driven in 
thermodynamic equilibrium state. The second part of investigation is intended to be studied in further publication. 
Here wc only discuss, at a qualitative level, how equilibrium state can be destroyed and point out possible mechanisms 
of transition to the turbulent state. 

As it was already mentioned, supposition about thermodynamic equilibrium state can be founded by considering 
vortex loops as a subsystem submerged into a thermostat with which they exchange energy. From the very fundamental 
physical principles it follows that a set of vortex loop should be driven in equilibrium state and, accordingly, to obey the 
Gibbs statistics. The underlined physical field plays a role of a thermostat in this case. In case of a Bose condensate, 
which we consider in given paper, the order parameter plays the role of such field. Vortex filaments in this study 

are described as crosses of the surfaces where both real and imaginary parts of the order parameter ■!/'(r,t) vanish. 
Other excitations of the order parameter (phonons, and rotons in case of real Hell) interact with an ensemble of 
vortex loops, driving the latter to thermodynamic equilibrium. In accordance with the principle of maximal entropy, 
the probability distribution function (functional) obeys the Gibbs formula. The given reasoning, however, bears a 
rather general character. It gives no details concerning dynamics of the system since they are hidden behind the 
temperature and chemical potential definitions. Therefore it does not answer the question: " how the system of 
vortices is driven in thermodynamic equilibrium? " Without knowing dynamical details we also can not to answer 
another question: " how the thermodynamic equilibrium state for a set of vortex loops is destroyed and the transfer 
to a strongly nonequilibrium (turbulent) state occurs? " We stress that we' are interested in analytical results, since 
there ia a number of numerical works (see [6], [16]), where the thermalization of vortices producing in equilibrium 
underlined field was demonstrated. 

It is well known however that Gibbs distributuion can be alternatively obtained in the frame of reduced models like 
kinetic equations or the Fokker-Planck equation (FPE) . These methods are not as general as a principle of maximal 
entropy, but instead they allows us to follow the details of how the Gibbs distribution is established. In this study, FPE 
was used to examine how chaotic vortex loops generated by the thermodynamic equilibrium Bose-Einstein condensates 
reach equilibrium. This study uses the Gross-Pitaevskii model [18] [19] (for details and notations see [23]). 

To achieve that goal we propose the following steps. First, we state a problem of stochastic dynamics for the order 
parameter V'(r,t). Dynamics of quantity ip{r,t) is supposed to satisfy the Gross-Pitaevskii equation [18] [19]. To impose 
the stochastic behavior of a BEG, the random stirring force /(x, t) was added to the right side of the Gross-Pitaevskii 
equation. The random force obeys the fluctuation-dissipation theorem. This guarantees that the Gibbs distribution 
for probability functional 7'({V'(r)}) = Af exp{—H {'ip{r)} /T) is followed. Here H {^{v)} is an energy of particular 
configuration {il){r)} of the order parameter, T -temperature, and A/" is a normalization factor. Then, starting from 
that dynamics, we derive an equation of motion for a vortex lines elements cosidering the latter as zeroes of the order 
parameter. The random force acting on the order parameter transforms into random velocities for elements of the 
vortex line. Then we derive the Fokker-Planck equation for probability distribution function functional 'P({s(^)}) 
in a vortex loop configuration space. Finally we show that the Fokker -Planck equation is satisfied by the Gibbs 
distribution exp(— _ff {s} /T), where iJ {s}- is functional of energy due to the vortex loops. Analyzing this result we 
discuss the possible reasons for destruction of thermodynamic equilibrium and suggest the possible mechanisms of 
transition to the turbulent state. 



I. LANGEVIN EQUATION FOR THE ORDER PARAMETER 

In this section we briefly describe the problem of the stochastic motion of Bose-condensate on the basis of Gross- 
Pitaevskii model [18] [19]. In the frame of reference where the normal velocity is absent the Gross-Pitaevskii equation 
reads 

= (1 - iA)(-|lvV -ni^ + Vo ivi' V). (1) 

Here A is (dimensionless) kinetic coefficient, n is chemical potential (energy per particle) and Vb is the potential 
describing an interaction of particles. The order parameter ^ is defined here so that 

= ^/p'se"^. (2) 
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The quantity ps is the superfluid density, (j) is (dimensionless) phase, the flow velocity v is defined as v = (h/m)W(f). 
Rescaling the time t t{h/m) and adding a random stirring force one can rewrite 1 in the following form: 



Here H {f/'j-is the Ginzburg-Landau free energy functional, having the following form 



The thermal noise obeys the following fluctuation-dissipation theorem 

(/(xi,ti)/*(x2,t2)) = 2kBTK 5(xi - X2) 8{h - t2) 



(3) 



(4) 



(5) 



The stochastic problem described above is the well studied problem, so called model A. It has a solution describing 
thermal equilibrium, where probability of some configuration of V-filed is proportional to exp(— {^'(r)} /T. It 
immediately follows from the correspondent Fokker-Planck equation [20]. 



II. LANGEVIN EQUATION FOR THE VORTEX FILAMENT 

In this section we derive equation of the stochastic motion of quantized vortex lines in Bose-condensate on the 
basis of the problem formulated in 3 and 5. Vortices arc crosses zeroes of the both real and imagine parts of the 
parameter, and the problem is just to describe the motion of zeroes of the (complex) function t/j varying in time and 
space. That problem, as well as other more general problems of that kind have been considered many times by many 
methods, starting from pioneering papers by Pitaevskii [21] and Fetter [22] (see also [23]). We however develop some 
new method, not too rigorous from mathematical point of view, but convenient for the main purpose of the paper. 

Vortices in quantum fluids are very slender tubes (except of the phase-transition region) and their dynamics is 
analogous to that of the strings. This implies that we have to aim our efforts at integrating out the radial degrees 
of freedom. It can be reached by the following procedure. As any motion of topological defects, the vortex line 
dynamics is determined by the underlying fleld theory. On the other hand, if one ignores for a while a presence of 
other excitations except of the vortices one can say that all motion of Bose-condensate is determined by a motion of 
(all elements) of (all) quantum vortices. Hence, there is a mutual correspondence and the order parameter tp{x,t) 
can bo considered as a fimctional of the whole vortex loop configuration ip{'x. \ s(^,t)). In particular, the temporal 
dependance of field ip is connected with the motion of lines and its rate of change (at some point x ) is expressed by 
the following chain rule. 



at J Ss{^',t) dt 

r 

Hereafter we will denote integrals along the individual lines and summation over all loops by the index F. Let us 
rewrite equation (3) in the following form: 



h? + l dt Stp* A? + V 

Let us multiply equation (7) by Sip* / 5s{S^q, t), where ^0 is some chosen point on the curve. Adding the result to the 
complex conjugate and integrating over whole space we have 



A-i d^p Sip* A + i d^* 6ip 



A'^ + 1 dt Ss{^o,t) A2-I-1 dt 5s{^o,t) 



J { sr <5s(eo,t) sr 5s(^o,t)y ^' 



A2 + l-'^ ' '5s{iQ,t) ' A2 + 1-' ^ ' '(5s(eo,i) 

The integral in the second line of (8) expresses a chain rule for functional derivative 5H{s)/ 5s(S^Q,t), where H{s) is 
the energy of Bose-condensate expressed via vortex line position. Consequently considering them to be very slender 
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tubes (which is justified when the radius of the curvature R is much larger of the core size tq), and neglecting the 
part of energy associated with the core, the quantity H{s) is just the kinetic energy of the superfluid flow created by 
vortices (see e.g. [24], [25]) 

r r' 

Here s'(^) is a tangent vector, k is the quantum of circulation equal to 27r?i/m. Calculation of the functional derivative 
dH{s)/ds{£,o,t) is straightforward and leads to result 

'^^^'^ =p,Ks'{^o)xB{^o). (10) 



Quantity B(^o) is the velocity of the vortex line element expressed with well known Biot-Savart law 

« /■ (s(eo)-s(0)xs'(0 

^^^'^--tJ |s(€o)-s(e)|3 (11) 

r' 

Thus, we expressed the integral in the second line of (8) purely in terms of variable s(^, t), describing the vortex line 
configuration. 

Now we will treat the integral in the Ihs of equation (8), which is the first line of this relation . Substituting Rel. 
(6) for temporal derivative ^ (and the complex conjugate expression) into integrand we have 

r' 

Integral (12) depends on the whole vortex loop configuration. Nevertheless this integral, (and other similar integrals 
including squared gradients of the order parameter) can be evaluated approximately in a general form, provided that 
the line is slightly curved (sec e.g. [23]). The relevant calculations, however, are rather cumbersome, and, to not block 
up the basic text, we have taken out them in the Appendix. As follows from the calculations given in the Appendix 
the term interesting for us is a just sum of (A. 5) and (A. 10), so the first line of equation ( 8) is transduced in 

^s(^o)xs'(6) + |?^s(eo). (13) 

Thus, we again expressed the considered term (first line of (8) in terms of variable s{^,t), describing the vortex line 
configuration. 

Let us now discuss the rest term of equation (8) including the random force /(x, t) (the third line). Consequently 
considering that all motion of Bose-condensate is connected to motion of line ( We remind, that we do not include 
acoustic efi'ects), we have to consider Langevin force /(x, f) as some secondary quantity stemming from random 
displacements of all elements of all loops. This implies that the small casual variations of the order parameter Sip are 
connected to the small casual displacements of lines with the help of the following chain rule: 

5V(x,t)= I gg|5s(^',i)dr. (14) 
r 

Furthermore, the casual small variations of Sip and small casual displacements of lines (for some interval of time 
St) are connected to random forces by an obvious fashion 

Ss{^,t) = f{^,t)St, (15) 



SiP{ic,t) = f{x,t)St. (16) 

The reader should distinguish random force /(x, t) acting on the order parameter and random velocity f (^, t), acting 
on the line on their arguments, besides f($,i) - is, certainly, vector. Substituting (15), ( 16) in (14), and reducing on 
St, we receive the following connection between random force f{x,t) and random velocity f{^,t): 
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= /2|r|f«'.«'. (17) 

r 

Further we should substitute this expression (and also conjugate expression) in the third line of equation (8), that 
gives 



r' 



(18) 



This expression coincides with expression (12) with an obvious replacement s (^o) — > f(^o,i)- Applying technique 
used in the Appendix we can directly write down the final result as the similar to expression ( 13) 

-f(eo,t)xs'(6) + -r^f(6,i). (19) 



Note that our consideration does not relate intensities of the random force /(x.t) acting on the order parameter and 
random velocity f(^,f), acting on the line, which will be done a bit later. Gathering together all transformed terms 
of equation (8), namely (10), (13) and (19), we have (after cancellation by 27r) 

^ s (^o) X s'(eo) + s (^o) - 



A2 + 1 ^^"^ ^^"^ A2 + 1 

= s'(Co) X B(Co) + ^J^^iM X s'(^o) + J^^M (20) 

Vector equation 20 can be resolved up to velocity S|| ^o) along the curve, which does not have physical meaning and 
can be removed by suitable parameterization of label variable ^. Resolving equation 20 we arrive at 

« (^o) = YTWB(eo) + ^^^^s'(Co) X B(^o) + f(^o, t). (21) 

Equation (21) describes the motion of vortex line in terms of line itself. It is the remarkable fact (not obvious in 
advance) that noise f(^0:i) acting on the line is also additive (docs not depend on line variables). 

The last effort which we have to do it is to ascertain both the statistic properties of noise f (^o , i) and its intensity. 
Equation ( 19) appeared as the result of transformation of the third line of equation (8) including random force /(x, t), 
so they are two equal vectors. Equating them we have 

^^f(C..)xs'(,,.A|^f(,.,.)= (22) 

= /'^^x.fA^/(x„.o^ + ^r(x..O: 



Taking a scalar productions of both parts of equation (22) on themselves taking at another time t2 and point ^2 we 
get 

47r2^^i^(f,,(a,ti)f.i(6,io) + f..(a,ti)fr,.(6,t2)) = (23) 



2 /■ /"jS j3 <5V'(xi,t2) ^V'*(x2,t2) , w*. ,^ 



Here f,,^ and f^^ are components of random velocities in ?7i,?72 directions in the normal (to the line) plain. We 
omit the rest terms in the rhs foreseeing them to drop out after averaging. Averaging relation (23) and using the 
fluctuation-dissipation theorem (5) we obtain 

(f,,(a,ii)f,.(6,i2)) = p^l^^i^-^ ^^l + K^X '^(^1-^2) 5{U - h)5r,,r,, . (24) 

Thus, starting from dynamics of Bosc-condcnsatc ( equation (3)) with fluctuation-dissipation theorem (5) we derive 
equation (21) describing motion of the vortex line in terms of line itself with the additive noise obeying (24). These 
relations complete a statement of the stochastic problem of quantized vortex dynamics under random stirring velocity 
stemming from the random force acting on the underlying field. In the next subsection we demonstrate that this 
problem has the equihbrium solution given by Gibbs distribution exp(— ff{s} /ksT),. Here -ff {s}- is the functional 
of energy due to the vortex loop (equation (9)) and T is the temperature of Bose-condensate. 
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A. The Fokker-Planck equation 



To show it we, first, derive the Fokker-Planck equation corresponding to Langevin type dynamics obeying (21) and 
(24). Let us introduce the following probability distribution functional 

n{<o},t) = {v^) = {H<o - ■ (25) 

Here (5 is a delta functional in space of the vortex loop configurations. Averaging is fulfilled over an ensemble of 
random forces. Rate of change of quantity V {{s{£^)} , t) is 

Here s^et is the deterministic part of velocity determined by ((21)). Omitting standard details of derivation (see for 
details e.g. [29]) we write the Fokker-Planck equation in a form 



dt J ^ (5s(0 

+ y d€ y de imm')) = o. (27) 

Here (f(^)f(^')) is the correlation function of random velocity applied to the line and satisfying (24). Our aim now 

is to prove equation (27 ) to have a solution in the form 'P({s(^)}) = Afexp(—H {s} /T), where A/" is a normalizing 

factor. We start the proof with the first term in integrand of (27). Using ( 10) and parameterization, that s ^o) is 
normal to the line, we write 

B(0 = ^/iO X (28) 
Let us rewrite the first term in integrand of (27) in the following form 

Performing the functional differentiation and using a tensor notation we rewrite (29) in the following form (we omit 
the coefficient before the integral): 

'«-<-^WA.r,.-£|^+ (30, 



/ 

+ J d^eM-H{s}/kBT)e"^'' (^^^^^^ 
+{l/kBT) J d^exp(-//{s}/A;Br)e"^^s^(0 



5 5H{s) 

(?)MCvO, 

SH{s) dH{s) 



The functional derivative 6s'i^{^) / Ssa{^) oc S^a, therefore all of terms vanish due to symmetry. Thus, the reversible 
term gives no contribution to the FP equation (27), one can say it is a divergent free term (sec e.g. [26]). It is clear, 
that the said above refers only to a case of thermal equilibriums, i.e. it is valid only for the Gibbs solution. 

The last term in the FP equation(27) can be transformed by use of fluctuation-dissipation relation for line (24) in 
the following way: 

J ^^'(f(Of(^'))^^= (31) 

Jd^l d^'^I^^!±^ S{^,-^,) d{t,-h)6r,,,r,, _i_exp(-iJ{s}/fcBT) (32) 

Observing (10) and (28) again one convinced himself that this term exactly compensate the second term in (27) as it 
should be. 

Thus, we have proved that the thermal equilibrium of Bose-condensate results in the thermal equilibrium of a set of 

vortex loops. In the next section we will review the basic steps of our derivation. On that ground wc will discuss how 
the thermal equilibrium can be destroyed and consider (qualitatively) a possible scenario of transition to a turbulent 
state. 
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III. ANALYSIS OF SOLUTION. POSSIBLE SCENARIO OF TRANSITION TO A TURBULENT STATE 



Analyzing our proof one can see that the following steps were crucial: 

1. The additive white noise /(x, <;) acting on field il){'x.,t) is transformed into the additive white noise f{£,,t) acting 
on vortex line position s{^,t). 

2. The intensity of noise /(^, t) expressed by (24) exactly compensates the dissipative flux of probability distribution 
functional in the Fokkcr-Planck equation (27). Furthermore, the compensation occurs in integrand, in other words 
the detailed balance between pumping and dissipation takes place. 

3. The non-dissipative part of equation of motion of the line elements gives no contribution to the flux of probability 
FP equation (27), it is divergent free (. [26]). 

Resuming the said above, it is possible to say, that by the basic reason for an establishment of the thermal 
equilibrium in set of the vortex loops is the detailed balance between primping and dissipation which occurs for all 
scales. It certainly agrees with the basic principles of the statistical physics (see, for example, [28]). It should be 
realized that the rest part of the equilibrium problem, namely evaluation of the partition function (and, accordingly, 
all physics) basing on the Gibbs solution 'P{{s{^)}) = A/'exp(— iJ {s} /T), is the far more difficult problem. The real 
success is achieved mainly in numerical simulations [6], [16]. As for analytical studies, we would like to point out 
a series of works by Williams with collaborators who uses 3D variant of the Kosterlitz-Thouless theory (for review 
and bibliography see [3]). We appeal to results of the cited works to continue our investigation. One of the most 
striking results obtained in these works is an observation that in equilibrium there are vortex loops of appreciably large 
sizes. That result can be hardly foreseen by naive estimation that probability of appearing of long loop is Boltzmann 
suppressed, ~ exp{—al/kT), where a is energy per unit of length. That implies, that, first, there is strong screening 
and, second, vortex loops have a polymer like structure with folding of line elements. Both phenomena reduce greatly 
the effective energy per unit of line length. Actually distribution of loops with respect to their lengths obey the law 
n{l) oc with some exponential factor (see [16]). 

We think, that the presence of the long vortex loops is responsible for possible destruction of thermodynamic 
equilibrium. Indeed, the small vortex loops (rings), of order a wave length of phonons, have quite definite energy 
and momentum. Degrees of freedom due to an extent of these objects are expressed weakly. In other words, the 
small vortex loops (ring) can be interpreted, as some kind of thermal excitations, which just renormalize a density 
of the normal component. Interaction with other quasi-particles (phonons and rotons) is reduced to usual collisions, 
that leads to a renormalization of usual kinetic coefficients. Quantitatively this renormalization is small, and can be 
appreciable only very close to the; lambda transition (see [3]). 

As opposed to the small loops, the long loops, first, have a very large amount of interior degrees of freedom. 
Secondly, the long weakly curved filament interact with the normal component under the laws applied to macroscopic 
rectilinear vortices (see, for example [1]). Thus, the long and small loops interact with a quantum fiuid in different 
manners. The situation here is quite similar to a 2D case, where an interaction with third sound arises only from 
the vortex pairs of the large sizes, while small vortex - antivortex pairs just renormalize the density of the superfluid 
component (see [31]). Because of a presence of long vortices, the detailed balance can be easily destroyed. Indeed, the 
strong enough large-scale perturbation (appeared e.g. due to nonuniform external stream) can violate the equilibrium 
state. Two scenarios, leading to the turbulent state are possible. First of them is connected to development of 
nonlinear Kelvin waves on the vortex line. Because of a nonlinear character of the equation of motion, the large-scale 
perturbations on initially smooth line interact, creating higher harmonics. They, in turn, produce higher harmonics, 
etc. In real space this process corresponds to the tangling of a long vortex loop and formation of parts with very 
high curvature (kinks or hairpins), which are subjected to the strong dissipation. Thus, in space of the Kelvin waves 
a nonequilibrium state is established. This state is characterized by the Kolmogorov flux of curvature along scale in 
the one-dimensional space describing the vortex lines (similar sc:(;narios are described in works [8], [9], [11]). 

Other scenario, leading to violation of the equilibrium state due to the presence of long loops, consists in the 
following. In the presence of the normal component the long vortexes are exposed to mutual friction, and, as a 
consequence to the Magnus force. As a resiilt, "ballooning" (in average) of the vortex loops is possible. Enlargement 
of an amount of the long loops displaces the law n{l) oc by slower dependence. Uncompensated long vortex 

loops are exposed to "redundant" reconnections, breaking down into smaller loops. This process can be interpreted, 
as the nonequilibrium Kolmogorov cascade in space of the loop sizes. The reader familiar with the theory of superfluid 
turbulence, easily finds out, that the second scenario coincides with the one, introduced by Feynman for nonequilibrium 
dynamics of the vortex tangle in the counterflowing Hell (see reviews [1], [7], and also original clause [27]). The only 
diflference is that we consider the given process on a background of thermal equilibrium, and the cascade arises from 
deviation from n{l) cx distribution. 

Certainly, both scenarios described above are not independent. For example, at the moment of reconnection there 
are strong deformations on smooth lines, generating the Kelvin waves. On the other hand, internal dynamics of the 
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vortex filaments influences frequency and details of reconnection. Thus, analyzing the received equilibrium solution, we 
have offered the scenarios, according to which, the transition to strongly nonequilibrium, turbulent state, is possible. 
We would like to stress, that realization of these scenarios is due to the presence of very long (actually macroscopic) 
vortex loops in the equilibrium solution. 

Let's discuss now, what sorts of the macroscopic perturbations are capable to destroy thermodynamic equilibrium. 
First of all, it would be nonuniformities appearing during a flow of supcrfluid in channels. Other case is the counterflow 
of the normal and superfluid components appearing in the heat flux. Superfluid turbulence can be also generated by 
intensive first and second sounds [32]. Furthermore, the vortices, which have arisen in fast transition throw lambda- 
point at the first moments are in nonequilibrium, and, probably, relax under the turbulent law (compare [4], [5]). 
There are other, more exotic, mechanisms of destroying of thermodynamic equilibrium, for example, with use of the 
ion current. 

In Introduction we have mentioned about the problem of initial vorticity, which was faced by Vinen in his description 
of the transient superfluid turbulence (see [14]). In order to make his equation to be self-consistent, he introduced the 
phenomenological, "generative" term, whose origin is not clear. Under the scenarios described above, the transition 
in turbulent state is just development of instability of long vortex loops with respect to the large-scale perturbations. 
Since instabilities develop very rapidly (usually exponentially), a presence of a plenty initial loops is not necessary. 
Thus, the phenomenological "generative" term in the Vinen equation is just some composition of a small amount 
of the long vortex loops and instabilities increments. Of course, the confirmation of that point of view requires the 
quantitative analysis of the problem of stability of vortex loops, in the flowing (counterflowing) superfluids, this is 
supposed to be studied in the further work. 



Exploring the Langevin dynamics of a chaotic set of quantized vortex filaments appearing in a weakly nonideal Bose 

gas, we have followed the mechanism of establishment of thermodynamic equilibrium. We have concluded, that the 
equilibrium state arises due to the additive random forces, acting on the order parameter, transforms into additive 
random velocities, acting on the lines in a proper manner, to satisfy the fluctuation-dissipation relation for the lines 
dynamics. It guarantees the detailed balance between random velocity and dissipation. Noting, further, that in the 
vortex tangle there present very long vortex loops, subjected to large-scale perturbations, we offer the qualitative 
scenarios of a destruction of the thermodynamic equilibrium and transition to the turbulent state. 
The author is grateful to E. Sonin and G. Williams for fruitful discussion. This work was partly supported by grants 
N 03-02-16179 from Russian Foundation for Basic research and grant N 2001-0618 from INTAS. 



A major contribution into integrals like (12) appears from a vicinity of the vortex filament. More precisely these 
integrals diverge, however a divergency has a weak logarithmic character and usually the radius of curvature (or 

interline space) is used as an upper cut-off. In practice this implies that we can replace functional x/'(x j s(^, t)) by the 
static vortex solution [23] ipv{^±) = ^i;(s(^ci)i) — x), where ^ci is the label of point of the fine closest to point x. It 
is clear that vector s{^cl, t) —x = x± lies on the plain normal to line and crossing point s{^cl,t). In accordance with 
the said above, the main contribution into integral comes from the points x remoted from the line at the distances 
much smaller than R, namely )s(^cj,<) — x| <C i?. The static solution tpy{x±) can be described as follows. The order 
parameter is V-u(x_l) = j4()xj_]) exp[i (/?], where A(|x^|) cx |x^| for |x_l| <C tq (core size) and A(|x^|) « y/pj for 
|x^| ^ tq. The phase (p is just a polar angle. It is also useful to write tp-uix^) in Cartesian coordinates ^71,7/2 lying 
on the normal (wrt line) plain and crossing points s(£^ci,t). Then ipyix^) = rji^ + -q^) cxp[i arctan 772//7i, ]. In 
vicinity of vortex it is approximately 'i/'i)(x_L) cx rji +i?]2- Thus, to evaluate integral (12) we replace ipi^i- \ s(^,t)) with 
the function ■i/'u(xj_) = V'u(s(^cZi i) — x) and the integrating over d^x with d^x_L(i^c;. Functional derivative 6ip* /Ss{^', t) 
should be evaluated by the following rule: 



IV. CONCLUSION 



APPENDIX: APPENDIX 



'^^(x|s(g,t)) 



Likewise, 



^V'(x|s(g,^)) 



V±^„(xi) ^(^0 
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Then the integration over d^cl and d^' removes delta- functions and the whole expression is taken in point ^o- Finally 

we have 



A-i 



A2 + 



J J (f^^{s (eo)ViV.(x±))ViV:(x±). (A.l) 



Here s (^g) is velocity of line at point of observation. The second term in integrand of the first line of equation (8) is 
just complex conjugate, which is 



A + i 



A2 + 



J J (f^^{s (^o)VxV:(x±))ViV.(x±). (A.2) 



The sum of (A.l) and (A.2) splits into parts with sum of integrands and their difference. The difference term can be 
rewritten with the help of famous formula of the vector algebra in form 

^ J rf'xx s (^o) X (VxC(x±) X (A.3) 



A2 

To evaluate (A.3) let us rewrite the static solution Vt)(xx) in local polar coordinates (|xj_| ,(p) Vector production 
Vx^*(x_l) X V_lV'i;(x_l) takes the form 

(Vx<(xx) X VxV.(xx)) = 2i^i^^^MMs'(eo), (A.4) 

x± 



And integral enterring (A.3) is evaluated as 

j d2x_L(V_LV«(x_L) X V_lVi;(x_l)) = 

= s'(eo) / 2i^^^i^^M^27r|xx|d|xx| = 27rzp,s'(eo). 
J |x_l| 

Thus the difference part of (A.l) and (A.2) represented by (A.3) takes the form 

2ttPs 



^2 + 1 « (^o) X «'(?o) (A-5) 
Let us now study contribution, which arises from the sum of (A.l) and ( A.2). It is 

^^2^ J cP^^ {(s (eo)VxV.(x±))VxV:(x±) + (s (?o)Vx^:(xx))Vx^.(xx)} . (A.6) 

The evaluation of (A.6) is easily fulfilled in Cartesian coordinates rji eta2 in x_l plain. Then the first term is 

(s (&)V.^„(x.))V.«(xx) = (s, ^) (Me. + ^e.) . (A.7) 

Here Si,S2 are components of the line velocity in 7717^2 directions, respectively, and ei,e2 are unit vectors in these 
directions. Note that expressions in first brackets of rhs of above relations are scalars whereas the second ones are 
vectors. Adding ( A.7) with complex conjugate, we obtain 



ei I 2 Si ^-|- S2 S2 I + 

V drji drji dr]2 drji drji dr]2 ) 

e2 {2 S2 si si 

V dr]2 dr]2 dr]2 drji dri2 drji J 



Combination 

dr]2 dr]i drji dri2 



(A.8) 
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(the second and the third terms in both brackets) disappears for static solution (xj_), because the complex conju- 
gation is equivalent to change r]i —r]2)- Due to the symmetry the §^§|f is just 5 | VxV'w(''^-l)I^ and integrating it 
gives the result: 



/ 



(fx± |V_lV^(x_l)|^ = 2'KPs In ^ = 2npsa, (A.9) 



where oo is some number (weakly depending on R). Resuming all said above we finally have that the term sought 
transforms into 

s (Co). (A.10) 

The sum of expressions (A. 5) and (A. 10) is the first line of the equation (8), written down in the terms of variable 
s(C,i), describing a vortex configuration. 
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